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ABSTRACT 


We reanalyze data of near-infrared background taken by Infrared Telescope 
in Space (IRTS) based on up-to-date observational results of zodiacal light, inte¬ 
grated star light and diffuse Galactic light. We confirm the existence of residual 
isotropic emission, which is slightly lower but almost the same as previously re¬ 
ported. At wavelengths longer than 2 /in i, the result is fairly consistent with 
the recent observation with AKARI. We also perform the same analysis using 
a different zodiacal light model by Wright and detected residual isotropic emis¬ 
sion that is slightly lower than that based on the original Kelsall model. Both 
models show the residual isotropic emission that is significantly brighter than the 
integrated light of galaxies. 

Subject headings: cosmology: observations — diffuse radiation — infrared 
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1. Introduction 

The extragalactic background light (EBL) has been observed over a wide range of 
wavelengths to examine the energy density of the universe. In particular, the near-infrared 
EBL has been thought to provide an important clue to our understanding of the early 


( Cambresv et al. 

2001; 

Goriian et al. 

2000; 

Wright and Reese 

2000 


Levenson et al. 

2007 

) and the InfraRed Telescope in Space (IRTS) ( 

Matsumoto et al. 

(2005 

), hereafter 


referred to as paper I) discovered that a significant fraction of the near-infrared isotropic 


emission cannot 


De exp 


(Tsumura et al. 


ained with known foreground emission. Recent AKARI observations 


2013dl ) also show a consistent result with COBE and IRTS at wavelengths 


longer than 2 /jrn. This excess background emission in the near-infrared sky is particularly 


interesting in light_of t 


sky flKashlinskv et al. 


re recent dis co very of large excess flue 


2005 


2007a 


bl; 


Matsnmoto et al. 


2011 


uation of the near-infrared 


Zemcov et al. 


2014). 


The results of paper I attracted wide interest due to the high accuracy enabled by the 
unique low resolution spectroscopy of IRTS, and its poin t sou rc e dete ction limit ( ~ 11 


mag) which was much deeper than COBE. However, as 


Mattila 


020061 ) pointed out, paper I 


did not take the contribution of the diffuse Galactic light (DGL) into account. Uncertainty 
of the zodiacal light (ZL) model has been also raised, since ZL is th e do minant f oregro und 


emission with a spectrum similar to the residual isotropic emission ( Dwek et al. 


2005). 


In response to these concerns, we decided to reanalyze the IRTS data using up-to-date 
observa tions of ZL, integrated star light (ISL), and DGL. While in paper I, we adopted the 


Kelsall et al. 


(119981 ) ZL model, in this work, we 


ZL model (the so called Wright model, 
two models. 


Wri g ht 


perform the same analysis using a different 


(I19981 U. and examine the difference between 


The overall outline of the paper is as follows. In section 1, we briefly present the IRTS 
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observation and the acquisition of raw data in section 2. In section 3, we estimate the 
contribution of foreground emission, ZL, ISL and DGL based on the latest observations. 
In section 4, we search for the residual isotropic emission for two ZL models based on the 
correlation with the sky brightness after subtracting the ISL and DGL. Finally, in section 
5, we discuss the astrophysical implications of the detected excess brightness. 


2. IRTS observations 


We now give a brief description of the IRTS mission and the data acquisition process. 
Details of IRTS mission can be found in paper I. 


IRTS was one of the mission experiments on the small space platform, Space Flyer Unit 
(SFU), that was launched on March 18, 1995. On a low-inclination near-earth orbit, IRTS 
continuously surveyed the sky avoiding both the sun and the eart’ 
lasted for about 30 days, during which 7% of the sky was surveyed i 


i. IRTS observations 


Murakami et ah 


1996). 


The Near InfraRed Spectrometer (NIRS) is one of the focal plane instruments of IRTS, 
and was optimized to obtain spectra of the diffuse 
of the flight performance of NIRS can be found in 


lackground ( 

IS 

oda et al. 

1994 

). Details 

Noda et al. 

( 

1996 

). NIRS covered a 


wavelength range from 1.4 /mi to 4.0 /mis with a spectral resolution of 0.13 /mi, providing 
24 independent wavelength bands. The beam size was 8 arcmin square and the detection 
limit for point sources was ~11 mag, and both are considerably better than those of COBE. 


The limiting magnitudes at wavelengths shorter anc 


the sky fluctuation and readout noise, respectively fjNoda et al 


longer than 2.5 /m i are dominated by 


( 19961 ). paper I). 


The NIRS detectors were of the charge integrating type and the ramp curves with one 
cycle of 65.54 sec were sent to the ground for all of the 24 wavelength bands. We retrieved 
the signal of the sky using 5 sec integration times (i.e. the signal difference in 5 sec in ramp 
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curve), during which no distinguishable stars and no cosmic ray hits were detected in any 
band. We obtained signals of the sky brightness after subtracting the dark current when 
the cold shutter was closed. During each 5 sec integration, the telescope axis moved about 
20 arcmin along a great circle, resulting in a trapezoidal beam pattern, 8 arcmin x 20 
arcmin in area. To make the contribution of stars and Galactic emission less effective, high 
Galactic latitude (b > 40°) data were extracted from the full data set for the background 
radiation analysis. The highest Galactic latitude was 58°, while the ecliptic latitude ranged 
from 12° to 71° in the selected sky. Complete spectra of the sky were secured at 1010 fields. 
Sky coverage of 1010 fields is ~ 60% of the surveyed area. 


3. Foreground emission 
3.1. Zodiacal light (ZL) 


The zodiacal light (ZL) is the emission component of the solar system which consists of 


scattered sun’ 


Kelsall et al. 


ight_and thermal emission by interplanetary dust. We adopted the model by 


( 19981 1 in paper I, which is a physical model constructed using the seasonal 


variation of the ZL observed wit li the Diffuse InfraRed Background Explorer (D1RBE) on 


COBE. At the same time, 


Wright 


( 19981 1 proposed a different physical model based on 


the so called ’’strong zodi principle” assuming no residual emission at 25 /jrn towards the 
ecliptic pole. In this paper, we use both models and examine the difference of residual 


i sotrop ic emission. As for the Wright model, we used the model revised by 


( 2000 ). 


Gorjian et al. 


We retrieve the brightness of the ZL in the DIRBE bands corresponding to the 1010 
IRTS fields at the epoch of the IRTS observations for both models, and construct the ZL 
model brightness for the IRTS bands. We obtain the model spectrum for the scattering 
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and thermal emission part, separately. For the scattering part, we simply normalize the 
spectral shape of the sun (ASTM G173-03 Reference System) [j] at the K band model 
brightness. For the thermal part, we extrapolate the M band model brightness to the 
shorter IRTS wavelengths, assuming a 300K blackbody in accordance with recent AKARI 


observations (ITsumura et al. 


2013b). The ZL model spectrum is obtained by summing the 


scattering part with the thermal emission part. Compared to the model used in paper I, 
the ZL component at wavelengths longer than 3 /mi is slightly brighter. In the L band, our 
adopted models render a ~ 6% brighter value than the original model, but still within the 
uncertainty of the models. The validity of the adopted ZL spectrum will be compared with 
observations in section 4. 


3.2. Integrated star light (ISL) 

The integrated star light (ISL) for stars fainter than the limiting magnitude composes 
a portion of the foreground emission. The best way to estimate the ISL is to sum the 
brightness of stars in the 2MASS catalog that fall within the beam, however, the uncertainty 


of the attitude determination and irregular beam pattern elongated al ong scan pat 


l makes 


this analysis difficult. While, in paper I, we applied the SKY model ([Cohen 


19971 ). and 


assume a simple cosec{b ) law to the model ISL at three selected fields, in this paper, we 

based on an improved model, that is, the TRILEGAL Galaxy model 


obtain the moc 

cl ISL 

(Girardi et al. 

2005) 


2005). We further calculate ISL for 12 equally-spaced fields along the scan 


path for the I, J, H, K^, L, and L^ bands. 

Since the fluctuation of the ISL is not negligible, we perform 50 Monte-Carlo simulations 


1 Available in electronic form at 


http://rredc.nrel.gov/solar/spectra/aml.5/astmgl73/astmgl73.html 
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for a 1 square degree field for all TR1LEGAL stars fainter than the IRTS limiting magnitude 
and assign their average values to be the TR1LEGAL 1SL. However, the overall uncertainty 
of the TR1LEGAL model is so large, that we must calibrate the TR1LEGAL ISL by 
comparing with the 2MASS ISL for the H and K bands. The 2MASS ISL is obtained by 
summing the ISL for the 2MASS stars fainter than the IRTS limiting magnitudes, and that 
for stars fainter than 2MASS limiting magnitudes (15.8 mag for the J band and 15.1 mag 


for the H band; 


Skrutskie et ah 


( 2006! )) which is estimated by adopting the TRILEGAL 


model. Contribution of the latter part to the 2MASS ISL is only a few percent. Fig.l shows 
the correlation between the TRILEGAL ISL and the 2MASS ISL for the 12 selected fields. 
Left and right panels indicate the results for the H and K bands, respectively. Horizontal 
errors in Fig.l represent the 1 a dispersion of the TRILEGAL ISL as a result of these 
Monte-Carlo simulations and indicate the expected fluctuation of TRILEGAL ISL. Fig.l 
clearly indicates that the 2MASS ISL is brighter than the TRILEGAL ISL both for the H 
and K bands, and the scatter of the 2MASS ISL is consistent with that expected from the 
TRILEGAL model. Based on this analysis, we finally obtain the model ISL by multiplying 
TRILEGAL ISL by 1.23 with 8 % error. The ISL for the 1010 fields observed with IRTS 
was obtained by interpolating the two neighboring model fields assuming a cosec(b) law. 
The ISL for 24 of the IRTS bands was estimated by interpolating the model ISL using the 
blackbody spectrum with the same limiting magnitudes as in paper I. 


Besides the model errors, we calculate the model ISL for variation of ± 0.5 mag in 
the limiting magnitude and assigned the difference to be peak-to-peak errors. This error 
is a little lower than the error due to uncertainty of the model ISL. The model ISL thus 
obtained renders a more reliable brightness and spatial distribution than that of paper I. 





3.3. Diffuse Galactic 


light (DGL 
as mentioned in 


The DG 


Mattila 


was not taken into account in paper I as a foreground emission, 


20061) . since no reliable observation of DGL was reported at that 


time and the contribution of DGL to the overall sky brightness was thought to be small. 
Recently new observations of the near-infrared DGL have been attained, and we attempt 


to estimate the DGL for the IRTS fields and bands. 


Tsumura et al. 


f 2013d ) obtained low 


resolution spectra of the diffuse sky with AKARI for the wavelengths ranging from 2 /mi 
to 5 /jin, detecting a clear correlation b etwe en the_near-infrared sky brightness and the 


far-infrared emission (100 /mi, 


Schlegel et ah 


( 19981) ). A PAH band at 3.3 g m was clearly 


etec 


;ed, however, the detection was limited to the Galactic plane at b < 15°. 


Arai et al. 


(120151) performed a simil ar co rrelation ana 


Spectrometer (LRS) (ITsumura et al. 


ysis with data from the Low Resolution 


2013al ). one of the instrume nts on the sounding ro cket 


experiment, CIBER (Cosmic Infrared Background ExpeRiment, 


Zerncov et al. 


(120131 )). 


They also detected a clear correlation with the far-infrared emission for the wavelength 
range from 0.95 /mi to 1.65 /mi at high Galactic latitudes. 

Fig.2 summarizes their results. The ratio of the DGL to the far-infrared (100 /mi) 
emission is shown in units of nWm“ 2 sr _1 / /iWm _2 sr _1 . Filled and open circles represent 
the result of GIBER and AKARI, respectively. The CIBER result indicates scattered star 
light by interstellar dust, and the AKARI result shows the thermal emission of the fine 
dust particles tran siently heated by a single UV photon. We fit the scattered part wit 


the ZDA04-BC03 (Brandt and Draine 


2012) model recommended by 


Arai et al. 


(12011 


as is shown with the dotted line. As for the thermal part, it is not cl ear t hat the AKAR 


result can be applied to the DGL at high Galactic latitudes, since 


Tsumura et al. 


(2013c) 


reported lower level of thermal DGL than expected from Fig.2 at higher Galactic latitudes. 
Therefore, we first assume that no thermal part exists at IRTS fields where Galactic 
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latitudes are higher than 40°. The thermal contribution, especially from the 3.3 //m PAH 
feature will be discussed again in section 4. Adopted values for IRTS bands are shown by 
open squares in Fig.2 for which an error of ±20% was applied. 


We used the far-infrared map compiled by 


Schle gel et al. 


(119981 1. and retrieved 100 


/irn brightness for a 12 aremin diameter FOV for the 1010 IRTS fields. Using the ratio of 
the DGL to the 100 /irn brightness (Fig.2), we are able to obtain DGL for 24 IRTS bands 
and for 1010 IRTS fields. 


4. Residual isotropic emission 

We attempt to obtain the residual isotropic emission with same procedure as paper 
I. First, we subtract the ISL and DGL from the observed sky brightness and make a 
correlation analysis with the model ZL brightness. The upper sets of data points in Fig.3 
show typical correlation diagrams at 1.8 /irn both for the Kelsall model (left panel) and the 
Wright model (right panel). For both models, we find an excellent linear correlation for all 
wavelength bands, and intersection of linear fit line at x = 0 provides the residual isotropic 
emission. The lower sets of data points in Fig.3 show the individual residual emission 
for the 1010 IRTS fields at 1.8 /irn after subtracting all foreground emission components, 
demonstrating that the residual emission is fairly isotropic. Since error levels are almost 
same for two models, it must be noted that there is no clear preference between these two 
models. 

Fig.4 shows the dependence of the residual emission at 1.8 /irn on the Galactic latitudes 
for the case of the Kelsall model. Data points indicate averaged values for 10 degrees along 
the Galactic longitude. Filled circles and open circles indicate the result of present work 
and that of paper I, respectively. The large scale structure observed in paper I disappeared, 
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and the data points of the present work show random scatter. This improvement is mainly 
due to our revision of the ISL model. 


Fig.5 shows slopes of the linear fit line both for the Kelsall and Wright model. In 
both models, slopes are a few percent larger than 1.0, which is consistent with paper I. We 
regard this as the deviation of the ZL spectrum from the solar spectrum which reflects the 
physical properties of interplanetary dust. The ZL spectrum can be obtained by multiplying 
the solar spectrum by these values, showing reddening of the ZL color at near-infrared 
wavelengths. The observed spectrum smoothly connects to the ZL spectrum from 0.75 to 
1.6 /iin observed with CIBER and is consistent with the_J band brightness of the Kelsall 


model based on the DIRBE data (see Fig. 9 in 


Tsumura et ah 


(201 


®). 


Fig.6 shows the residual isotropic emission obtained for both the Kelsall and Wright 
model, with two different sets of error bars. The inner and outer error bars indicate the 
random and total error, respectively. Aside from the random error, the systematic error 
makes the spectrum change in the same direction. Random errors include: fitting errors 
from correlation analysis, ISL errors due to limiting magnitudes, calibration errors, and 
DGL errors. Systematic errors are due to model errors of ZL and ISL. As paper I, errors in 
the ZL model are e stimated by interpolating the uncertainties in the original model (Table 


7 in 


Kelsall et al. 


( 19981 1). ZL model error is the dominant source of error and amounts to 


~ 80 % of total error. Table 1 indicates numerical values for the residual emission and their 
error for two models. 


Compared with residual emission in paper I, which is based on the Kelsall model, the 
peak brightness of the residual emission of the present work is ~ 11 nWm _2 sr _1 lower. The 
residual emission at wavelengths longer than 2 /irn is almost the same as that of paper I. 
A flat spectrum for the three shortest wavelength bands is a characteristic feature found 
in this analysis. The residual emission obtained by adopting the Wright model for the ZL 
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provides ~ 7 nWm _2 sr _1 lower peak value than that obtained by adopting the Kelsall 
model but the spectral shape is very similar to each other. Fig.6 implies that there exists 
excess near-infrared isotropic residual emission independent of the choice of the two ZL 
models used here. 


Additionally, we make the same analysis including the thermal part of the DGL. In 
this case, residual emission at 3.28 /mi PAH feature is estimated to be ~20 % (~ 2.6 
nWm _2 sr _1 ) lower than that of Fig.6, causing a sharp absorbing feature in the spectrum of 
the residual emission. We also find that the sky brightness at 3.28 /mi after subtracting 
ZL and ISL shows no correlation with far-infrared sky brightness. These results favor the 
non-existence of the PAH feature at high Galactic latitudes, however, this is not conclusive 
given that random noise at 3.28 /mi is so large. In any case, the contribution of the thermal 
part of DGL is almost negligible and does not significantly change the final result. 


Fig. 7 shows the breakdown of the emission components whose spectra are obtained 
as the average brightness at high ecliptic latitudes (/? > 70°) and high Galactic latitudes 
(b > 45°) for the case of the Kelsall model. The ZL Spectrum indicated is obtained by 
applying the result of correlation analysis (Fig.5). The residual isotropic emission amounts 
to ~ 1/4 of the ZL, and this is comparable with the seasonal variation of the ZL. 


Fig.8 summarizes the observations of the residual isotropic emission for which the 


Kelsall model is appli ed for the ZL. Typical two observations flCa/mbresv et al 


Levenson et ah 


2001 


2007) are shown for COBE data. Solic 


integrated light of ga 


deep Galaxy counts flKeenan et ah 


axies (I LG ) by 


Totani and Yoshii 


ine indicates the model of the 


(120001 ). which is consistent with 


2 0101 ). The brightness of this work is a little 


than that in paper I, but still consistent with COBE and AKARI (Ts umura et aL- 


ower 


2013d). 


It must be emphasized that these three satellite observations render consistent residual 
emission within error, although the beam size and limiting magnitudes are completely 
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different. In the optical region, two new observations are plotted. Matsuoka et ah (2 0111) 
r e-an aly zed the Pioneer 10/11 data, where the ZL is negligible near the Jupiter orbit, while 


Mattila et ah 


(20111) attempted to detect the EBL using shadowing effect of dark cloud 


at high Galactic latitude. Both results show that the observed background brightness is 
consi ste nt with the known foreground emission. On the other hand, HST observations by 


Bernstein 


(120071) show fairly bright residual emission at a similar level as this work. 


Fig\9 indica 


Bernstein 


;es e xcess b right ness o ver ILG. Since some galaxies are already removed in 


(12007 b and 


Tsumura et ah 


( 2013d ). we subtracted 10 and 50 % of the ILG (solid 


line in Fig. 8) from their data in Fig.8, respectively, based on ILG magnitude relationship 


by 


Keenan et ah 


(120101) . We find improved consistency between IRTS and AKARI data. 


As already mentioned in paper I, the blue Rayleigh-Jeans like spectrum is clearly seen at 
wavelengths longer than 1.6 //m. 


5. Discussion 


It has been suspected that the residual isotropic emission is a part of t 


theispectrum of the residual isotropic emission is similar to that of the ZL ijPwek et_ah 


re ZL since 


20051 ). However, it is not so easy to construct the new ZL model which includes the 
residual isotropic emission, since the residual emission is comparable to the seasonal 
variation of the ZL. One possible approach is to add the new dust component which has a 
heliocentric and spherically symmetric distribution. It is, however, difficult to maintain the 
spherically symmetric distribution of interplanetary dust of inner solar system due to the 
perturbation by planets. Maintaining the supply of dust against Pointing-Robertson drag 
is also difficult. These considerations imply that the new dust component must be beyond 
the orbit of Jupiter. If we attribute the residual isotropic emission to the scattered sunlight 
by interplanetary dust, the brightness extrapolated to visible wavelengths amounts to ~ 
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100 nWm 2 sr 1 , assuming a solar spectrum. Th e exist ence of such a 


0/11 observation s (IHanner et al 


Matsuoka et al. 


component, however, 


1974). A recent reanalysis 


(2011) confirms that the residual isotropic 


is not detected with Pioneer 
of the Pioneer 10/11 data by 
emission near the orbit of Jupiter is less than 10 nWm _2 sr _1 at 0.44 and 0.64 pm. These 
results indicate that the new dust component does not exist, otherwise the new component 
dust has peculiar optical properties. 

The light of the first stars at the re-ionization epoch is another possible emission source 
which has been thought to be an important clue to delineate the star formation history 
of the universe. However, recent theoretical works based on high redshift galaxies predict 


much l ess co ntribution of the first stars to the near-infrared background flCoorav et al. 


2012b; 


Yue et ah 


2013 


It has been thought that the spatial fluctuation of the sky directly provides the 


characteristic feature of the EBL, since the fluctuation of the ZL is so low (Pyo et al. 


2012). Large fluctuations at angles larger than 100 arcsec that cannot be 


explained with the 


2005, 

2007a. 

o) and AKA 

II at 2.4, 3.6 and 4.2 pm ( 

Matsumoto et al. 

201 

-1). 

Jecent Spitzer 

) observations 

(Kashlinskv et al. 

2012; 

Coorav et al. 

2012a 

) and AKARI ( 

Seo et al. 

2015 


confirm that flat fluctuation spectra extend to degree scales. The spectrum of fluctuation is 


blue Rayleigh-Jeans like (IMatsumoto et al 


2011) and clear spatial correlations are found 


for both AKAR1 and Spitzer wavelengths. 


Zemcov et ah 


(2014) recently report the result of sounding rocket observations with 


C1BER. Th ey detect large fluctuation at 1.1 and 1.6 pm. In Fig.9 we over-plot the result of 


Fig.2 in 


Zemcov et ah 


(120141 1 as large filled squares. Unit is shown in right ordinate with 


same dynamic range as excess brightness. The spectrum of the fluctuations is similar to 
the excess brightness, suggesting that the excess fluctuation and brightness may be of same 
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origin. A rat io of the absolute brightness to the fluctuation, //A/, is ~ 30. Furthermore, 


Zemcov et ah 


(20.14) find excellent spatial correlation among 1.1, 1.6 and Spitzer 3.6 /mi. 


Since no clear signature of redshifted Lyman a is detected, they conclude that the first 
stars are not the source of ex cess f luctuation. They also exclude direct-collapsed black holes 


(DCBHs) model flYue et al. 112013m . s ince the detecte c 


is quite different from the prediction. 


Zemcov et ah 


fluctuation is too high and the color 


12014) claim that the large portion 


of the observed fluctuation can be explained with intra halo light (IHL) and estimate the 
surface brightness of IHL to be comparable to that of 1LG. It must be emphasized that the 
quoted EBL (ILG plus IHL) is still a few times lower than that of the residual emission 
obtained in this work. 


At present, there is no definite emission source which can explain the observed excess 
brightness. More observational and theoretical work is needed to delineate the emission 
source of excess brightness and fluctuation. 


6. Summary 

The reanalysis of the IRTS data was performed to obtain improved and more reliable 
measurement of the near-infrared residual isotropic emission. We revised the estimation of 
the integrated star light due to faint unresolved stars and the thermal emission part of the 
zodiacal light spectrum. We take the diffuse Galactic light into account based on the recent 
observations with AKARI and GIBER. Besides the Kelsall model for the zodiacal light used 
in paper I, another model, the Wright model, was examined, too. 

The result for the Kelsall model shows peak value of 60 nWm -2 sr -1 which is ~ 

11 nWm _2 sr _1 lower than that of paper I. This is still considerably brighter than the 
integrated light of galaxies (ILG). The result for the Wright model shows peak value of 53 
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nWm -2 sr -1 which is slightly fainter than that of the Kelsall model. Both models render 
significant residual isotropic emission that cannot be explained with known foreground 
emission sources. 
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Integrated star light, TRILIGAL (nW.m 2 .sr ’) 


Integrated star light, TRILIGAL (nW.m 2 sr’ 1 ) 


Fig. 1.— The correlation between the TRILEGAL ISL (the integrated light for stars fainter 
than the IRTS limiting magnitude based on the TRILEGAL Galaxy model) and the 2MASS 
ISL (the integrated light for 2MASS stars fainter than IRTS limiting magnitude and for 
stars fainter than 2MASS limiting magnitude based on the TRILEGAL Galaxy model) for 
the 12 selected IRTS fields. Left and right panels show the case for the H and K band, 
respectively. Horizontal error bars represent the standard deviation of the TRILEGAL ISL 
which was obtained using Monte Carlo simulations. Straight lines show the best fit for linear 
correlation, and the dotted lines indicate the adopted ±lcr error. 
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Wavelength (|im) 


Fig. 2.— The ratio of the DGL to the far-infrared (100 /jrn) emission, shown in units 


of nWm 2 sr 1 (DG 


GIBER ( Arai et al. 


)/ |iW m 2 sr 1 (FI R). Filled and open cir cles represent the result of 


2015) and AKARI ( Tsumura et 


shows the model (ZDA04-BC03 in 


Brandt and Draine 


2013c), respectively. Dotted line 


(120121 )) recommended by 


Arai et ah 


(20151). and open squares indicate adopted ratio for IRTS bands. 
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Zodiacal light, model brightness (nW.m 2 .sr ') Zodiacal light, model brightness (nW.m 2 .sr ') 


Fig. 3.— Upper set of data points shows the correlation diagram between the surface bright¬ 
ness after subtracting ISL and DGL from observed sky brightness at 1.8 /mi, and model ZL 
brightness. Left and right panels indicate the case for the Kelsall and Wright model, respec¬ 
tively. Solid lines are best fit lines for linear correlation. Lower set of data points represents 
individual residual emission after subtracting all foreground emission from observed sky 
brightness in which solid lines show residual emission obtained by linear correlation analysis. 
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Fig. 4. — Dependence of residual emission at 1.8 /ini on the Galactic latitudes for the case 
of the Kelsall model. Filled circles and open circles represent the result of this work and 
that of paper I, respectively. Data are takes by averaging individual residuals for 10 degrees 
along Galactic longitude. 
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Fig. 5.— Wavelength dependence of the slopes of the linear fit lines in the upper part of 
Fig.3. Slopes represent ratios to the solar spectrum. Filled and open circles represent the 
case for the Kelsall model and the Wright model, respectively. 
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Fig. 6.— Spectra of residual isotropic emission after subtracting ZL, DGL and ISL are shown 
for the Kelsall (filled circles) and the Wright model (open squares), respectively. For clarity, 
the residual emission of the Wright model is shifted to slightly longer wavelength. Two sets 
of error bars are plotted for each data point. The inner bars represent random errors, while 
the outer bars indicate the total error including systematic errors. 
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Fig. 7.— A breakdown of the sky brightness at high ecliptic latitudes ((3 > 70°) and high 
Galactic latitude (b > 45°) for the case of the Kelsall model. Filled circles, bars, filled 
squares, open circles, and open diamonds indicate the observed sky brightness, zodiacal 
light (ZL), residual isotropic emission, integrated star light (ISL) and diffuse Galactic light 
(DGL), respectively. 
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Fig. 8.— Summary of observations of the residual isotropic emission for the case of the Kel- 


sall model. Symbols represent fo llowing data: filled circles 
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IRTS, this work 


(12013. 
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, fillec 
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(1200111 ). open 


(2007)), filled diamonds 


(201111). and open diamonds (dark cloud, 


Mattila et al. 


(20111)). The solid line shows the integrated light of galaxies (ILG) (Totani and Yoshii 


2000 ). 
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Fig. 9.— Summary of the excess emission over integrated light of galaxies (ILG). Sym- 


bols represent t h e follow ing data: filled circles (I RTS. this work), h 


Tsumura et al. 
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Table 1: Surface brightness of residual isotropic emission and its errors in units of nWm -2 sr -1 


Wavelength 

[/H 

Kelsall 


Wright 

Residual 

Total Error 

Residual Total Error 

3.98 

16.0 

4.0 

15.6 

4.1 

3.88 

15.8 

3.6 

15.4 

3.7 

3.78 

12.9 

3.3 

12.8 

3.4 

3.68 

10.6 

3.7 

10.9 

3.9 

3.58 

14.0 

3.0 

13.9 

3.1 

3.48 

12.7 

3.1 

12.6 

3.3 

3.38 

12.9 

3.0 

12.7 

3.2 

3.28 

11.6 

3.1 

11.6 

3.3 

3.17 

15.1 

3.0 

14.9 

3.2 

3.07 

18.5 

3.0 

18.2 

3.3 

2.98 

17.1 

3.2 

16.6 

3.4 

2.88 

19.0 

3.5 

18.2 

3.7 

2.54 

22.3 

4.2 

20.9 

4.4 

2.44 

20.5 

4.6 

19.2 

4.9 

2.34 

23.6 

4.9 

21.7 

5.1 

2.24 

29.2 

5.3 

26.8 

5.6 

2.14 

35.4 

5.9 

32.2 

6.2 

2.03 

38.0 

6.8 

34.4 

7.2 

1.93 

40.1 

7.8 

36.3 

8.2 

1.83 

42.4 

8.7 

42.8 

9.0 

1.73 

53.1 

10.1 

47.3 

10.6 

1.63 

58.3 

11.8 

51.4 

12.2 

1.53 

59.9 

12.7 

52.8 

13.3 

1.43 

58.1 

13.1 

51.2 

13.7 






